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Extragalactic Planetary Nebulae (PNe) are not only useful as distance signposts or as 
> ' tracers of the dark-matter content of their host galaxies, but constitute also good indicators 

of the main properties of their parent stellar populations. Yet, so far, the properties of PNe in 
the optical regions of galaxies where stellar population gradients can be more extreme have 
remained largely unexplored, mainly because the detection of PNe with narrow-band imaging 
or slit-less spectroscopy is considerably hampered by the presence of a strong stellar back- 
ground. Integral-field spectroscopy (IFS) can overcome this limitation, and here we present a 
study of the PN population in the nearby compact elliptical M32. Using SAURON data taken 
with just two 10-minutes-long pointings we have doubled the number of known PNe within 
the effective radius of M32, detecting PNe five times fainter than previously found in narrow- 
band images that collected nearly the same number of photons. We have carefully assessed the 
incompleteness limit of our survey, and accounting for it across the entire range of luminosity 
values spanned by our detected PNe, we could conclude despite having at our disposal only 
15 sources that the central PNe population of M32 is consistent with the generally adopted 
shape for the PNe Luminosity Function and its typical normalization observed in early-type 
galaxies. Furthermore, owing to the proximity of M32 and to ultraviolet images taken with the 
Hubble Space Telescope, we could identify the most likely candidates for the central star of a 
subset of our detected PNe and conclude that these stars are affected by substantial amounts 
, of circumstellar dust extinction, a finding that could reconcile the intriguing discrepancy pre- 

viously reported in M32 between the model predictions and the observations for the later 
stages of stellar evolution. Considering the modest time investment on a 4m-class telescope 
that delivered these results, this study illustrates the potential of future IFS investigations for 
the central PNe population of early-type galaxies, either with existing SAURON data for many 
more, albeit more distant, objects, or from campaigns that will use the future generations of 
integral-field spectrographs that will be mounted on 8m-class telescopes, such as MUSE on 
the Very Large Telescope. 
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1 INTRODUCTION 

In the field of extra-galactic astronomy, Planetary Nebulae 
(PNe) are perhaps regarded mostly either as useful ind i cators 
for the distance of their galactic hosts dCiardulloetal.1119891 : 
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Ijacobv, Ciardullo, & Fordlll990Ujacobv et al.lll992l) , thanks to the 
almost universal - though not fully understood - shape of their 
luminosity function (PNLF, generally in the [O III]A5007 emis- 
sion), or as tracers of the gravitational potential in the ou tskirts 
of galaxies dRomanowskv et al J [20031; iDouglas et alj[2007h . Yet, 
as summarised in the review of ICiardullol [20061) . PNe can also 
be used as probes of their parent stellar population. For the clos- 
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est or brightest PNe, the detection of critical but weak diagnos- 
tic emission lines such as [O III]A4363 allows to directly measure 
the temperature of the ionised-gas and thus the metallicity of PNe, 
which in turn makes it possible to co nstrain the chemical enrich- 
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ment history of their host galaxy (e.g., Richer, Stasinka. & McCalll 
Il999l : IJacobv & Ciardullolfl999l : iDopita et al.lll997h . In more dis- 
tant galaxies, it is still possible to place useful constraints on 
the Oxygen abundance of PNe by using brighter lines such a 
[O ra]A5007, H/3, [N n]AA6548, 6584, Ha or [S n]AA6713, 6731, 
or to investigate the star-formation history of a galax y by study- 
ing the shape and normalisation of the PNLF (e.g. , iMarigo et al.l 
|2004 ISchonberner et alj|2007k iMendez et al"1l2008b . In fact, un- 
derstanding the origin of the PNLF is a puzzle that, once solved, 
promises to reveal new clues on the late stag es of stellar evolution 
and on the formation of PNe themselv es (e.g. JCiardullo et alj2005l ; 
iBuzzoni. Arnaboldi. & Corradill2006l) . 

In this context it is important to note that most known extra- 
galactic PNe have been found in the outskirts of their hosts, 
whereas the PNe population of the optical regions of galaxies re- 
mains largely unexplored. Most PNe studies have indeed been car- 
ried out through narrow-band imaging (e.g.JCiardullo etal.lll989h 
or slit-less spectroscopy (e.g.. |Pouglas et al.ll2007T) . where the de- 
tection of PNe is considerably hampered by the presence of a stel- 
lar background. Yet, it is in central regions of galaxies where stel- 
lar population gradients can be more extreme and where the stel- 
lar ages, metallicities and element abundances can be the most di- 
verse between different galaxies. Integral-field spectroscopy can 
overcome the previous instrumental limitations since it allows for 
the careful modelling of a galaxy's integrated stellar spectrum, and 
this paper aims to demonstrate the potential of IFS for studying the 
properties on PNe in the optical regions of galaxies using SAURON 
observations for the compact elliptical M32 (NGC221), our closest 
early-type galaxy. 

This work is organized as follows. In ij2] we briefly review the 
acquisition and reduction of the SAURON data for M32. In ij3]we 
detail how we optimised the extraction of the nebular emission in 
the SAURON data and how we proceeded to identify and measure 
the [O III]A5007 flux of the PNe in M32. In ij4]we present our main 
results, assessing in particular whether our data are consistent with 
the generally adopted shape for PNLF and the value for its normali- 
sation that is most commonly observed in early-type galaxies. In ij5] 
we discuss these findings by linking them to the known properties 
of the stellar population of M32, and in particular to the apparent 
dearth of po s t-asy mptotic giant branch stars that was reported by 
iBrown et alj d2000l) . Finally, in ij6] we draw our conclusions and 
consider some future prospects for more IFS studies of the PNe 
populations in the central regions of galaxies. 




Figure 1. Map of the equivalent width (EW) of the estimated [O III]A5007 
flux from narrow-band imaging based on the SAURON data. The flux in 
the [O III] region is measured through a 2.5A-wide spectral window cen- 
tred on the expected position of [O III] given the systemic velocity of 
V S ys = —197km s _1 for M32, whereas the flux of the stellar contin- 
uum is computed using the mean spectral energy density across the en- 
tire spectrum times the same 2.5A wavelength interval. The difference be- 
tween these two flux measurements gives an estimate of the [O III] flux, 
which once divided by the mean spectral energy density in the stellar con- 
tinuum leads to the mapped EW values. Isolated patches of positive EW 
values could originate from the unresolved [O III] emission of PNe, and in 
fact the brightest and most circular sources correspond to the 7 PNe that 
were a lready identified in the optical regions of M32 by Ciard ullo et alj 
11989L shown and labelled in black). The red circles delimit galactic re- 
gions that are most likely devoid of any nebular emission and where the 
SAURON spectra were added up to obtain spectra of very high quality. In 
turn, these high-quality spectra were used to derive the optimal templates 
that we adopted to describe the stellar continuum in our data (see text for 
more details). 



2 OBSERVATIONS AND DATA REDUCTIONS 

M32 (NGC221) was one of the special objects that were ob- 
served over the cour se of the SAURON representative survey 
dde Zeeuw et alj[2002h . and more specifically during the last run 
of that observing campaign at the 4m William Herschel Telescope 
and foll owing the installation of a new volume phase-holographic 
grating (Emsellem et al. 20041). The central regions of M32 were 
observed with two offset pointings during 600s for each expo- 
sure, while using the low-resolution mode of SAURON that gives 
a field of view of 33'.' x 44'.' fully sampled by "94 x 0"94 
squar e lenses (for more details on the instrument see Bac on et"al] 
l200ll) . The data from each pointing were reduced similarly to the 



data obtained for the objects of the main SAUR ON sample (see 
lEmsellem et al.1120041 ; iFalcon-Barroso et al-lfeoOol) . and the result- 
ing datacubes were merged and resampled in 0'.'8 x f/.'8 spatial ele- 
ments each corresponding to spectra with a final spectral resolution 
of4.2A (FWHM). 

The present SAURON data for M32 were alr eady used to ex- 
tract t he stellar kinematics that was modeled by Cappell ari et al.l 
(2006) and the only difference with the data used in that work and 
other papers of the SAURON proj ect is that here we did no t per- 
form any Voronoi spatial binning dCappellari & C opin 2003). This 
was done to avoid swamping the signal of the weaker PNe against 
an increased stellar background, and to allow for a more consistent 
analysis of the flux distribution from each of the unresolved PN that 
we may detect. 
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Figure 2. Map of the values of the A/rN ratio of the [O III]A5007 emis- 
sion measured in the SAURON spectra. The A/rN values are shown in a 
logarithmic scale and dark-blue bins identify regions where A/rN < 3 
and the [O III] line is not formally detected. Several of the weak positive 
patches in Fig.[T]correspond here to circular sources with gently declining 
A/rN profiles as would be expected from atmospheric scattering of the 
[O III] emission of an unresolved PN. This is not the case for the rectan- 
gular patch approximately 20" to south and east of the centre, which is a 
spurious measurement where in fact no [O III]A4959 emission is detected 
along with the larger [O III]A5007 line. 



3 DATA ANALYSIS 

3.1 Emission-line Measurements 

In order to identify PNe in M32 and measure their flux in the 
[O III]A5007 line we first need to separate as accurately as possible 
the stellar and nebular contribution to each of the SAURON spectra 
that sample the cen tral regions of this galaxy. For this purpose we 
used the method of ISarzi et al.l d2006l hereafter Paper V following 
the notation of the SAURON project) whereby a set of stellar tem- 
plates and Gaussian emission lines are fitted simultaneously to the 
spectr^E while following also the approach of lSarzi et al.l d20 ldh to 
further improve the match to the stellar continuum and ensure that 
the ionised-gas emission is extracted from the subtraction of a phys- 
ically motivated stellar model. This is achieved by using, instead of 
standard template libraries based on stellar spectra or single-age 
stellar population models, a more appropriate set of empirical tem- 
plates that are constructed while matching a number of high-quality 
SAURON spectra extracted from regions in the target galaxy where 
no ionised-gas emission is found. 

Fig. Q] shows how we identified such emission-line free re- 
gions in M32 using an [O III] narrow-band image made from 
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Figure 3. Map of the flux of the [Om]A5007 emission from M32, in 
10 _16 ergs _1 cm — 2 arcsec — 2 and in a logarithmic scale, together with 
the photometric measurements of both firmly and marginally detected PNe. 
Each of the PNe is labelled and delineated by red contours (in solid and dot- 
ted lines for solid and marginal detections, respectively) corresponding to 
the best-fitting Gaussian model to their observed flux distribution, which is 
used to compute their total [O III] flux -F5007 and corresponding detection 
limit (see text). The inner contour around each PNe shows the half-peak 
flux level of the Gaussian model, thus corresponding to a circle a FWHM 
in diameter, whereas the second contour shows the region containing 90% 
of the flux of the model around each PN, or of both the Gaussians models 
in the case of blended PNe sources. 



the SAURON data themselves (see the caption of Fig. [T] for 
details). Each of the spectra extracted from the circular aper- 
tures shown in Fig. 1 was fi tted with the pixel-fitting code of 
ICappellari & Emselleml d2004l) . over the full wavelength range of 
the SAURON data and u s ing th e entire MILES template library of 
Sanchez-Blazqu ez et alj ( 120061) . Both the quality of such aperture 
spectra and of our fit to them is very high, with values for the ratio 
between the median level of the spectra (S) and the average level 
of the fit residuals (residual noise, rN) ranging from 100 to over 
400, with S/rN ~ 300 on average. The weights assigned to the 
MILES stellar templates during each of these fits were then used to 
combine the MILES spectra into optimal templates that, owing to 
the excellent quality of our fit, can be in practice regarded as M32 
spectra deconvolved from the line-of-sight kinematical broadening. 
Fig. [T] also illustrates how our emission-free apertures are evenly 
spread over the SAURON field of view, which ensures that their 
corresponding optimal temp lates can accoun t for the presence of 
stellar population gradients dRose et alj|2005l) when they are used 
to match each of the single spectra sampling the central regions of 
M32. 

Fig.Q]also locates the PNe found in the narrow-band survey of 



Ciardullo et al. ( 1989), and reveals already the presence of weaker 
and isolated peaks of 5007A flux that may also originate from the 
unresolved [O III] emission of PNe. If this is the case, then the 
nebular spectra observed where the flux of an unresolved PN is 
scattered by the atmospheric point-spread function (PSF) should 
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Figure 4. Radial profiles of the flux of the [OlII]A5007 line (in 10 _16 crgs _1 cm~ 2 arcsec~ 2 ) of the PNe sources shown in Fig. [3] along with their 
corresponding best-fitting Gaussian model for the flux distribution (red lines, models are resampled in the SAURON 0'.'8 X 0'.'8 bins). Radial distances are 
computed from the center of the Gaussian models and the data points are color coded according to the value of the A/rN ratio. Colors change from blue to 
red for increasing values of A/rN till they saturate to white for A/rN > 6. Green corresponds to the detection threshold of A/rN = 3, above which the 
data points are also plotted together with error bars for the [O III] flux. The black lines show the flux of the Gaussian model that would fit the data if these were 
rescaled until the minimum value of A/rN within the FWHM (for a radius of 0'.'96) of the model reaches the threshold of A/rN = 3. When integrated, the 
best-fitting and rescaled Gaussian models yield the total flux -F5007 of each PN source and its corresponding detection limit, respectively. The model profiles 
for PNe 14 and 15 are shown with dashed lines to indicate that these are doubtful sources according the previous definition of the detection limit. The radial 
profiles extend out to a radius r = 3<rpsp, where crpsF was determined from the fit of the brightest source PN 1. PNe 3 and 7, 6 and 11,9 and 10 have been 
modeled simultaneously. 



be characterised by large values of the [O IIl]/H/3 ratio, virtually no 
[N I]AA5197, 5200 emission, and by very narrow lines that would 
not be resolved in the SAURON spectra. The standard fitting strategy 
of Paper V, whereby the [O III] doublet is fitted first and the H/3 flux 
is subsequently extracted holding the H/3 line profile to that of the 
[O III] lines, is therefore appropriate for isolating the emission of 
PNe, and here we only further imposed a fixed width on the [O III] 
lines, for an observed velocity dispersion of 108 km s _1 that cor- 
responds to the spectral resolution of the SAURON data (§ 2). 

Fig. [2] shows the map of the ratio between the best-fitting 
Gaussian amplitude A for the [O III]A5007 line and the noise level 
rN in the fit residuals. In Paper V we concluded that a minimum 
value of A/rN — 4 was required to detect [O III] emission, al- 
though this was while letting the width of the lines free to vary. 
Since instead here we have fixed the [O III] line width, our fits in- 
volve one less free parameter and we can measure unbiased [O III] 
fluxes also down to a A/rN — 3 threshold, similarly to the case of 
H/3 line measurements of Paper V when the [O III] line profile was 
imposed on the H/3 line. With such a detection limit, Fig.[2]confirms 
the presence of several isol ated weak patches of detected [O III] 
emission beside the PNe of ICiardullo et al. I lfl989h . and shows the 
absence of a diffuse ionised-gas component that is otherwise com- 
monly found in early-type galaxies (Paper V). 



3.2 PNe Detection and Flux Measurements 

To measure the flux of any known or possible PNe in the SAURON 
field-of-view, and at the same time establish whether a weak patch 
of [O III] emission is indeed consistent with the flux distribution ex- 
pected from an unresolved PN, we fitted each of the [O III] sources 
in Fig. [2] with a circular Gaussian function that is meant to repre- 
sent the shape of the PSF of our observations. The full-width at half 
maximum (FWHM) of such a PSF was first measured by matching 
the m ost cleanly detected o f our sources (corresponding to source 
21 o f lCiardulloetaill 19891) , and subsequently held constant when 
measuring the flux of other known or candidate PNe. The Gaussian 
model was resampled in SAURON C/.'S x 0'.'8 bins before compar- 
ing it to the data, and for close PN sources care was taken to de- 
blend their fluxes by fitting simultaneously two or more Gaussian 
functions. Based on these fits, we deemed a given source of [O III] 
emission a detected PN if A/rN > 3 within at least the FWHM of 
the Gaussian model. The model itself was then used to measure the 
[O III] flux of the PN, by integrating the entire Gaussian function 
out to radii where otherwise the [O III] emission generally would 
not be detected or could come from other nearby PNe. 

According to the previous detection criterion we find 13 PNe 
in the central regions of M32 observed by SAURON, with two ad- 
ditional sources that barely miss the detection threshold due to just 



PNe in the central regions ofM32 5 




4900 5000 5100 5200 
X (Angstrong) 




4900 5000 5100 5200 
X (Angstrong) 




4900 5000 5100 5200 
X (Angstrong) 



4900 5000 5100 5200 
X (Angstrong) 



160 PN7 




4900 5000 5100 5200 
X (Angstrong) 



4900 5000 5100 5200 
X (Angstrong) 




4900 5000 5100 5200 
X (Angstrong) 



4900 5000 5100 5200 
X (Angstrong) 



4900 5000 5100 5200 
X (Angstrong) 



4900 5000 5100 5200 
X (Angstrong) 



4900 5000 5100 5200 
X (Angstrong) 





150 




140 




130 






a 

tU 


iao 


-o 




X 


110 








100 




90 




140 




130 








120 


B 




0) 

•o 


110 


■A 




j3 


100 




90 




m 




4900 5000 5100 5200 
X (Angstrong) 



4900 5000 5100 5200 
X (Angstrong) 




4900 5000 5100 5200 
X (Angstrong) 



4900 5000 5100 5200 
X (Angstrong) 



Figure 5. SAURON aperture spectra for the PNe sources shown in Fig. [3]along with our best fit for the PNe emission and the background stellar continnum. 
Such spectra were extracted within a radius r = 3 (tpsf from the center of the best-fitting Gaussian models plotted in Fig. [4] and correspond to the sum of all 
the SAURON spectra observed in the bins shown in each of the panels of that figure. In each of the panels shown here the red line shows the sum of the stellar 
and nebular model (green and blue lines), with flux densities in 10 _16 ergs — 1 cm — 2 arcsec — 2 A _1 , whereas the small points show the fit residuals around a 
zero level that has been rescaled for clarity. The values of the A/rN ratio for the [O III]A5007 line is printed at the bottom of each panel, and show that the 
[O III] emission is detected in each spectrum, including those extracted around the doubtful PN sources 14 and 15. 



one value of A/rN being less then 3 within the FWHM of the 
models, which turned out to be l'.'94 for a ctpsf = 0'.'82. Each of 
these sources is identified and labelled in Fig. [3] which maps the 
[O III] flux observed in M32 and shows with red contours the best 
fitting Gaussian models to each of the detected PNe, including the 
questionable sources with dotted lines. The fit to the [O III] flux 
distribution of these sources can be better appreciated in Fig. [4] 
where the radial profile for the [O III] flux of each PN is plotted 
together with its best-fitting Gaussian model. In Fig.|4]the [O III] 
data points are colour-coded according to the value of the A/rN 
ratio measured in each bin to illustrate that down to our [O III] de- 
tection threshold of A/rN = 3 the [O III] fluxes follow very well 
the PSF model, which provides an independent check of our A/rN 
detection criterion for the [O III] emission. 

Fig.|4]also shows how for each source we determine the mini- 
mum flux that we could have measured at the particular position 
of the PN. In fact, we can expect such detection limits to vary 
across the galaxy, and more specifically to increase towards the 
center where only the brightest PNe lead to [O III] lines of suffi- 
cient strength to stick out above the higher level of statistical noise 
that comes with a brighter stellar continuum. Following our detec- 
tion criterion for PNe, we set the detection limit at the location of 
each our PN sources as the flux of the Gaussian model (shown with 
a black line in Fig. [4} that would match the [O III] fluxes if these, 
and the corresponding values of the [O III]A5007 line amplitude A, 
would be rescaled until the minimum A/rN value inside a FWHM 
reaches the threshold of A/rN = 3. As expected, for the doubt- 
ful sources the detection limits exceed the observed flux of the best 



model, since the minimum A/rN value is already less than 3 and 
therefore the data have to be scaled up, instead of down (although 
by just ~7% in both cases). The significance of the detection of our 
PNe sources determines the order in which the PNe are labelled in 
Figs. [3] [4] and subsequent figures. 

To further illustrate the degree to which all the PNe shown 
in Fig. [3] are detected, for each source we added all the SAURON 
spectra within a radius r = 3 upsf from the center of the best- 
fitting Gaussian models (corresponding to all the bins plotted in 
the panels of Fig. [4}, and fitted the resulting spectra with our stan- 
dard procedure. Fig. [5] shows such aperture spectra and our best 
fit to both the nebular emission from each PN and to the stellar 
continuum observed along the line-of-sight towards them. We note 
that the A/rN values for the [O III] emission measured in these 
spectra always indicate a detection, even within a 3 opsf aperture 
that is 6.5 times wider than the FWHM region where we decided 
to assess the detection of PNe and within which we might have ex- 
pected the emission from the weakest of our PNe to be lost against a 
larger stellar background. That this is never the case, however, sug- 
gests that the two questionable sources might as well be regarded as 
marginal detections, which is why we will keep considering them 
in the remainder of the paper. 

To conclude this section in Tab. Q] we list, for both firmly and 
barely detected PNe in M32, the position relative to the center of 
the galaxy, the total -F5007 flux of the [O III] emission with its corre- 
sponding detection limit, the velocity of each source, and, when the 
H/3 line was also detected, the average value of the [O IIl]/H/3 ratio. 
The latter two measurements are based on fits to spectra similar to 
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Table 1. Basic Properties of the PNe in the Optical Regions of M32 



PNe Luminosity Function of M32 



ID 


x-off 


y-off 


-F5007 


^5007, lim 


VpN 


[O IIIJ/Hp 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


1 


-4.48 


-15.50 


136.7 


7.6 


-214.2 


11.5 


2 


-5.16 


6.18 


101.5 


6.1 


-240.9 


13.8 


3 


-8.13 


1.91 


96.3 


6.9 


-190.2 


6.3 


4 


1.81 


-23.60 


46.5 


4.4 


-216.1 


7.8 


5 


-0.74 


3.86 


152.3 


14.5 


-260.0 


13.3 


6 


-4.75 


-5.83 


132.9 


15.9 


-240.4 


16.8 


7 


-9.75 


1.84 


76.2 


10.2 


-225.6 


7.7 


8 


5.99 


-3.29 


29.9 


6.4 


-188.9 




9 


9.78 


-6.45 


13.1 


5.6 


-188.1 




10 


12.18 


-4.38 


9.4 


6.0 


-147.5 




11 


-4.20 


-2.09 


16.7 


11.3 


-155.6 




12 


-11.33 


-9.69 


15.1 


11.2 


-133.6 




13 


8.92 


-20.59 


9.1 


8.2 


-189.1 




14 


4.08 


13.19 


7.9 


8.5 


-231.6 




15 


1.15 


-1.62 


52.6 


56.6 


-148.1 





Notes: (1) PN ID. (2)-(3) R.A. and Declination offset position, in arcsec- 
onds, from the center of M32. (4)-(5) Total [OlII]A5007 flux and detec- 
tion limit, in 10 _16 ergs — 'cm' 2 . (6) Velocity in km s — 1 , as measured 
in spectra extracted within a FWHM wide aperture around the center of the 
best-fitting Gaussian models shown in Fig. [3] (7) [O III]/H/3 in the same 
spectra, when H/3 is det ected. As shown i n Fig. [TJ our sources 1 to 7 were 
already detected by Ciardullo et al. ( 1989) who, in the order, labelled them 
with the ID of 21, 26, 24, 20, 27, 23 and 25. 

those presented in Fig. [5] but now extracted within a FWHM-wide 
aperture (rather than within a radius r = 3 ctpsf) in order to max- 
imise the emission-line signal and better isolate the kinematics of 
PNe that are close to each other in projection. 



4 RESULTS 

The analysis described in the previous section has delivered the 
solid detection of 13 PNe within the optical regions of M32, with 
an additional 2 sources where the observed [O III] flux distribution 
is only marginally consistent with the emission from an unresolved 
PN. Given the systemic velocity V sys = — 197km s _1 of M32 
and its average stellar velocity dispersion wi thin one effective ra- 
dius cr e — 60 km s _1 (Cappellari et al. 200rJ), the velocities listed 
in Tab. Q] show that the observed PNe very likely belong to M32, 
rather than to the disk or hal o of M3 1 where PNe move on average 
at a speed of -400 km s _1 dMerrett et"ai]|2006l) . In fact, the case 
for membership in M32 holds even at a local level, when the PNe 
velocity VW is compared with the mean stellar velocity V* and ve- 
locity dispersion <r* measured along the line-of-sight towards their 
location. Using FWHM-wide aperture spectra to better separate the 
velocity Vpn of blended or close PNe and larger 3 opsF-wide aper- 
tures to extract robust V* and ov measurements also in the outskirts 
of M32, we found that only in 3 cases | Vpn —V*\ > ct*, even less 
often than what is normally expected for 15 sources. Tab. [T] also 
lists, when we could measure it, large values for the [O IIl]/H/3 ra- 
tio, thus showing that these are indeed high-ionisation sources most 
likely consistent with PN s pectra. Finally, we note that for the PNe 
that were also detected by ICiardullo et al.l dl989t) our flux values 
agree fairly well, within a few percent, with the values reported in 
that work. For instance, for the brightest PN in our sample (source 
5) we m easure a F5007 flux tha t is only 8% fainter than what was 
found bv lCiardulloe7aT] ( ll989l source 27 in their Tab. 7). 
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Figure 6. Luminosity function of the PNe sources observed with SAURON 
in the optical regions of M32, including marginal detections (dashed gray 
histogram) and along wit h the theoretical form of the PNe luminosity func- 
tion (red solid line, from Ciardullo et al. 1989). The absolute magnitude of 
each PNe is shown at the top of the figure with filled circles that are plotted 
at an arbitrary constant ordinate and that are connected to open diamonds 
showing the magnitude corresponding to the detection limit of each source 
(with marginal sources shown in grey). The dashed line shows the PNe 
luminosity function multiplied by the median values of the completeness 
function across the entire SAURON field, whereas the dotted lines indicate 
the range by which the expected number of detected PNe would vary de- 
pending on their exact position within the 0'.'8 X 0'.'8 SAURON bins (see text 
and Fig.|7]for details). By integrating such a completeness-corrected lumi- 
nosity function we obtain the total number of PNe that we would expect to 
detect within the region mapped by SAURON, and by matching this value 
with the actual number of observed PNe, 15 including marginal detections, 
we obtain our best estimate for the normalisation of the PNe luminosity 
function of M32. At first sight, the observed PNLF would look inconsis- 
tent with the theoretical prediction even accounting for incompleteness, in 
particular at its brightest end. A Kolgomorov-Smimov test and the simula- 
tions shown in Fig. [^demonstrate however that such a discrepancy can be 
expected with only 15 objects. 



Fig. [6] presents the luminosity function of the PNe found in 
the optical regions of M32 by our SAURON observations, which es- 
sentially map this galaxy within its effective radius Re = 30" . To 
constr uct such a PNLF we followed the definition of lCiardullo et all 
(1989) to compute the apparent V-band magnitude 7715007 = 
—2.5 log -F5007 — 13.74 of our sources, and derived the values for 
the absolute magnitude M5007 that are shown in Fig. [6] by adopt- 
ing a distance modulus of 24.49 magnitudes. This corresponds to a 
distance of 791 kpc, and it is the same value based on the surface- 
brightn ess fluctuation measurem ents of lTonry et al.ld200ll) that was 
used by Cappella ri et al. (2006). 

Fig. |6] also show s with a red lin e the expected form for the 
PNLF introduced by ICiardullo et all who combined the simple 
iHenize & Westerlundl 1 1963b model for a population of slowly fad- 
ing and expanding PN envelopes with a sharp exponential cutoff 
at the high-mass end. The shape of our observed PNLF is incon- 
sistent with such a theoretical curve, but this is hardly surprising 
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Figure 7. Reconstructed optical images of M32, in logarithmic flux scale, showing with grey bins the regions where, from left to right, PNe of increasing 
absolute magnitude M5007 and decreasing brightness would not be detected. To check whether in a given bin of the SAURON field a PN of absolute magnitude 
M50Q7 would be detected, we generated at that position a Gaussian model for the [O III] flux of a PN of that brightness in M32, deriving also the corresponding 
spectral density values for the amplitude of the [O III] line. Using the values for the level of the noise rN in the residuals of our spectral fits we then computed 
the values of the A/rN ratio around the PN position and, using the criterion introduced in § 13.21 simply checked whether A/rN > 3 within a FWHM from 
the centre of the Gaussian model, which in these particular simulations correspond to the center of the SAURON bins. 



considering that at the faint end our PN number counts are af- 
fected by incompleteness. To a first approximation the onset of 
this bias can be appreciated by noticing in Fig. [6] how the values 
for the absolute magnitude corresponding to the detection limits of 
each of our sources (shown with open diamonds) pile up from a 
M5007 ~ —1.5 and brighter , corresponding to an apparent mag- 
nitude limit of m5007 ~ 23. In general, while looking at bigger 
galaxies than M32 with a correspondingly larger number of PNe or 
at very close objects such as the large and small Magellanic clouds, 
other studies use the PNe that are found within the completeness 
limit in order to determine both the best normalisation for the the- 
oretical PNLF and to test whether this could in fact be considered 
the parent distribution for the data. In the case of M32, however, 
there are simply not be enough PNe within the completeness limit 
(say for A/5007 < — 2, taking a conservative guess from Fig. [6] to 
confidently carry out such a measurement and test. We therefore 
decided to first understand how our observa tional biases would af- 
fect the theoretical form of ICiardullo et al.l for the PNLF, to then 
test whether all our data, over the entire M5007 range, could have 
been drawn from the corresponding completeness-corrected model 
PNLF, while also determining the best normalisation for it. 

Fig. [7] presents simulations designed to show how we derive 
the completeness function of our observations, that is, the prob- 
ability as a function of absolute magnitude M5007 that a PN of 
that brightness could be detected across the entire field-of-view of 
the SAURON observations for M32. Since at any given luminos- 
ity the number of expected PNe scales with the total number of 
stars, such a probability is just the fraction of the total stellar flux 



encompassed by our field-of-view that is observed within the area 
over which PNe of a given absolute magnitude M5007 could be 
detected. To then decide whether a PN of given M5007 could be 
detected at any given position, we generated a Gaussian PN model 
of total flux corresponding to the apparent magnitude 7715007 and, 
using the values of rN obtained from our spectral fits across the 
entire field-of-view, simply checked whether A/rN > 3 within a 
FWHM of the Gaussian model. As Fig. [TJ illustrates, starting from 
an absolute magnitude of M5007 ~ —2.5 (msoo7 ~ 22) it becomes 
more and more difficult to detect PNe near the central region of the 
galaxy, with the area of non-detectability (shown with grey bins) 
quickly expanding outward for PNe fainter than M5007 ~ —1.0 
(fi5007 ~ 23.5). More specifically, we are 100% complete down 
to M5007 ~ —3, whereas we lose about 5% of PNe with an ab- 
solute magnitude of M5007 ~ —2.5 and nearly 26% of those with 
M5007 ~ —1.0. Fig.|7]shows simulations for just eight M5007 val- 
ues and while placing the Gaussian PN models at the center of each 
SAURON Of.' 8 x Of! 8 bin, but to derive the completeness function and 
the uncertainties associated with it we adopted a more refined grid 
of M5007 values and accounted for the impact of randomly placing 
the Gaussian models within the SAURON bins. 

When the theoretical PNLF of Ciardul lo et alj is multiplied by 
the completeness function that we have just derived, we finally ob- 
tain the completeness-corrected prediction for the expected num- 
ber of PNe in M32 that is shown by the red dashed curve in Fig. [6] 
with the dotted red lines tracing the scatter in this correction due 
to the exact PNe position in the SAURON bins. A Kolgomorov- 
Smirnov test reveals that this corrected PNLF can be regarded as 
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the parent distribution for our data, since there is a 82% proba- 
bility that all of our 15 PNe could have been drawn from it, and 
a 52% probability when considering only the 13 secure PNe. By 
integrating the completeness-corrected luminosity function we ob- 
tain the total number of PNe that we would expect to detect, and 
by matching this value with the actual number of observed PNe, 
15 including marginal detections, we obtain our best estimate for 
the normalisation of the PNe luminosity function of M32 within 
the region mapped by our SAURON observations. Such a normal- 
isation is generally expressed in terms of the luminosity-specific 
number density of PNe, through the ratio a = Np-^ e / Lboi be- 
tween the total expected number of PNe iVpNc in a stellar popu- 
lation and the total bolometric luminosity Lboi of the latter. Typ- 
ically, the total num ber of PNe is estimated by integrating the 
normalised PNLF of ICiardullo et alJ from the brightest observed 
absolute magnitude M*,5 qq7 = —4.47 to 8 magnitu des fainter, 
which is the limit where iHenize & Westerlundl d 19631) locate the 
faintest PNe. Yet, since the faint end of the PNLF has not been 
well constrained by observations, in some instances a is provided 
only within the completeness limit and such an extrapolation is 
avoided. Most often in the literature this limit lies within 2.5 mag- 
nitudes of M*,5oo7 and the luminosity-specific number density pa- 
rameter is therefore indicated as 02.5. Adopting our normalisa- 
tion there should be 76 PNe in total within the central region of 
M32 mapped by SAURON and between 7 and 8 PNe (7.6) with 
•W5007 — Af*,5007 ^ 2.5. Given that our observations cover M32 
within essentially its effective radius R e from I-band images, we 
are here encompassing nearl y half of the total I-ba nd luminosity, 
i.e. Li = 2.3 x 1O 8 M (see lcappellari et al.l2006l for R e and to- 
tal luminosity I-band measurement s). Adopting the reddening cor- 
rected B — V = 0.88 c olour of Buzzoni. Arnaboldi. & Corradil 
(2006) for M32, and using iBruzual & Chariot! d2003h models for a 
Salpeter IMF and Solar metallicity, we then obtain a B — I = 2.01 
and thus a B-band luminosity of Lb = 1.3 x 1O 8 M . Finally, 
following lBuzzoni. Arnaboldi. & Corradil J2006l) also for the bolo- 
metric correction to the B-band luminosity we arrive at a bolo- 
metric luminosity Lboi = 3.4 x 10 8 M Q , which means that for 
the optical regions of M32, within R e , the luminosity-specific PN 
number density is a — 2.2 x 10 — 7 Lq , and a factor ten less for 
0:2.5. Considering a 33% statistical error associated with our nor- 
malisation based on 15 points (for an intrinsic number of 25 PNe 
within our detection magnitude ran ge), these values are remarkably 
consis tent with thefindings of bof h lBuzzoni. Arnaboldi. & Corradil 
J2006l) and lCiardullo et al.l dl989h who give a = 1.70 x lO^L" 1 
and 0:2.5 = 2.23 x 10 _8 Lg 1 , respectively. 



To conclude, we have detected 15 PNe in the optical regions 
of M32 that were mapped by our SAURON observations, and thanks 
to the integral-field nature of our data and by carefully accounting 
for the incompleteness of our observation we have been able not 
only to confirm kinematically that these sources belong to M32 but 
also that their observed luminosity function and their total num- 
ber are consistent with the g enerally adopted shape of the PNLF 
(that of ICiardullo et al.ll 19891) and the typical values for the PNe 
number density a in early-type gal axies (a few 10 -7 per L , 
iBuzzoni. Arnaboldi. & Corr adi 2006). Considering that our 600s- 
long SAURON observations at the 4m William Herchel Telescope 
collected nearly as many p hotons as the 3600s- and 1800s-long 
narrow-band observations of lCiardullo et al.l d 19891) at the 0.9m and 
2.1m Kitt Peak telescopes, respectively, the fact that we more than 
doubled the number of known PNe in the central regions of M32 
while reaching 5 times fainter [O III] PNe fluxes (see TabQJ illus- 



trates well the power of integral-field spectroscopy to study PNe in 
the optical regions of galaxies. 



5 DISCUSSION 

It is interesting to comment on our PN results in light of what is 
known about the central stellar population of M32 and the current 
expl anations for both t he shape and normalisation of the PNLF. 

iRose et all d2005h have shown that the luminosity-weighted 
mean stellar population of M32 at lR e is older by ~3 Gyr and 
more metal-poor by about —0.25 dex in [Fe/H] than the central 
value of ~4 Gyr and [Fe/H]~ 0.0, which are trends that conspire 
to produce the flat color profiles previously reported in this object 
dPeletier|[l993l : I Lauer et ai1 Fl998). Episodes of star formation may 
enhance for a relative short period of time the bright-end of the 
PNe luminosity function of a galaxy, since massive main-sequence 
stars tend to produce also m assive and bright cen tral PNe stars 
(see e.g., the case of M33 in ICiardullo et1iuT 2004). For instance, 
lMarigoetal] ( l2004l) shows that stars with initial mass between 2.5 
and 3.5M evolve to central PN stars of about 0.70 - O.75M that 
can power nebular fluxes corresponding to the typical PNLF cut-off 
magnitude M*,5007 of -4.47. Since massive stars are short-lived, 
however, such a PNLF enhancement would be pretty weak already 
3 Gyr after star formation ceased, and it would no lo nger apply to 
the br ightest PNe but to objects with M5007 ~ —2.5 jMarigo et al.l 
120041) . It is thus not surprising, in particular given the small num- 
ber of PNe at our disposal, that we do not see such a signature in 
our observed PNLF. If anything, there seems to be a lack of PNe 
around that intermediate M5007 value. 

On the other hand, and like in other early-type galaxies with 
no evidence of recent (less than 1-Gyr-old) star formation, the pres- 
ence in M32 of bright PNe at a time when the progeny of massive 
stars have long disappeared is a puzzle that has not yet been fully 
solved. In fact, the need for a massive central star in the brightest 
PNe, with mass above ~ 0.70Mq, is a result that still holds in the 



studies of Schon berner et al.l d2007l) and iMendez et al.l d2008l) that 
supersede the initial investigation o f lMarigo etal.ld2004l) . One pos- 
sible way to form the required high-mass PN cores, suggested by 
ICiardullo etai]d2005l) . involves the binary coalescence during their 
hydrogen burning phase of ~ IMq stars, a process that has been 
associated with the formation of blue-stragglers and which could be 
facilitated in early-type galaxies by the abundance of ~ IMq stars 
in old stellar systems. If such a second path to the formation of 
PNe exists, then we may also expect to find a PNLF characterised 
by a peak at high M5007 above an otherwi se monothonically in- 
creasing PNLF, corresponding to a standard IHenize & Westerlund 
population of PNe powered by only old central PN stars of rela- 
tively low mass. In light of this scenario, it would be tempting to 
see a bimodality in our observed PNLF (Fig. |6), which happens 
to show two peaks around A/5007 ~-3.5 and -1 that correspond 
to parent stellar masses near 2 and 1M , respectively. However, 
due to small number statistics we cannot read too much in the ob- 
served shape of the M32 PNLF, as it is instructive to see through 
simple simulations such as those presented in Fig. [8] which con- 
firm the formal consistency indicated by a simple KS-test between 
our data and the completeness-corrected ICiardullo et al.l form of 
the PNLF. Similarly, it would be difficult with our sample size to 
claim a significant difference in the kinematic behaviour of bright 
and f aint PNe (as in NGC 4697, ISambhus. Gerhard & Mended 
2006), which could occur for instance if the coalesced cen- 
tral stars of the PNe can form in globular clusters, as could 
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Figure 8. Comparison between the spatial distribution and the luminosity function of the PNe detected in the central region of M32 (top and lower left panels, 
respectively) and three particular synthetic realisations of the same quantities (central and right panels). More specifically, the top panels show only the spatial 
distribution of PNe that were or would be detected in the optical regions of M32 by displaying values for the A/rN ratio above the A/rN > 3 detection 
threshold, and which correspond (adopting the observed values of the residual-noise level rN) to the flux distribution of either our be st-fitting Gaussian 
models to the PNe of M32 ( S!3.2t or to similar Gaussian flux profiles for randomly generated PNe. Starting from an intrinsic Ciardullo et al. shape of the PNLF 
(red curves in lower panels) and adopting a normalisation leading to match the number of observed PNe once the PNLF is corrected for incompletess (red 
dashed lines, see fj4j, the synthetic PN fields were generated by considering at any particular position in SAURON field of view the probability of having a 
PN of a given luminosity. Such a probability function simply corresponds to the total intrinsic PNLF rescaled by the fraction of stellar light that is observed 
in the SAURON spatial bin that is being considered. The [O III] flux of the simulated PNe was then "observed" by obtaining maps for the A/rN ratio and by 
applying the same detection criteria described in £|3.2| These simulations illustrate that, given the low number of PNe expected in M32, not much can be read 
in the apparent bimodality of the PNLF that we observe nor in the fact that the PNe of M32 seems to all be on the same side of the galaxy. In fact, these three 
particular simulations yielded the same total number of detected PNe than we observe (15), but generally this is not the case and the PNLF of such synthetic 
fields appear to differ even more from the observed PNLF than shown by the chosen examples. 



be the case for X-ray binaries (Sarazin, Ir win & Bregmanl koOO; 
IWhite, Sarazin & Kulkarnil2002h and where indeed blue-stragglers 
are often found. Alte rnatively, and consistent with the hypothesis of 
iMarigo et all d2004l) . a different kinematics for the brightest PNe 
may be a sign of recent accretio n of a younger and smaller galaxy 
(Mam on. Dekel, & Sto"eiu1l2005h . Yet, integral-field data allow for 
a direct comparison with the stellar kinematics as a function of 
M5007, such as shown for M32 in Fig. [9] that would be interest- 
ing to apply to more massive systems hosting a larger number of 
PNe. 

From a theoretical perspective, our best normalisaton of the 
PNLF in the optical regions of M32, with a value for the specific 
PN number density a in line with that of other early-type galax- 
ies, also poses a problem. Indeed, as long as the time-scale for PN 
visibility depends mostly on the lifetime of their post-AGB cores 
(pAGB), then since cores of lower mass spend more time in their 
high-temperature regime a should increase with the stellar pop- 
ulation age, whereas in fact quite th e opposite is observed (see 
iBuzzoni. Arnaboldi. & Corradill2006l for a comprehensive review 



of this discrepancy). A possible solution to this puzzle is to as- 
sume that a considerable fraction of the stellar population of early- 
type galaxies ends up in the extreme hot end of the helium-burning 
horizontal-branch (HB), and that when such stars of very small 
hydrogen-envelope mass subsequently leave the HB they not only 
skip entirely their AGB phase (thus the name of AGB-manque) but 
they also fail to produce a PN. As the AGB-manque evolution can 
effectively transfer some fraction of the pAGB energy budget to the 
integrated UV flux of the galaxy, this scenario could also explain 
the anti correlation between th e strength of the UV-upturn (as de- 
fined by B urnstein et al.l 19881) and the luminosity-specific PN num- 
ber density of early-type galaxies, which respectively increase and 
dec rease with bo t h gala xy mass and metallicity. In fact, it is in M32 
that lBrown et al.l (2000) provided some support to this picture, find- 
ing with near-UV Hubble Space Telescope (HST) observations that 
the relative faint UV flux of M32 originates from hot HB stars and 
their AGB-manque progeny rather than from their brighter post- 
AGB counterparts. Yet, given the relatively long UV-bright lifetime 
of stars leaving the hot end of the HB (~ 10 6 — 10 7 yrs compared 
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Figure 9. Velocity Vpn of the PNe in M32 relative to the average stellar 
velocity V* along the line-of-sight towards the PNe position and normalised 
by the stellar velocity dispersion cr* in the same direction, sorted by their 
absolute magnitude Msoo7- PNe sources are labelled and colour-coded as 
in Fig [6] The dotted horizontal lines for |Vpn — V*|/cr + = 1 delineate 
the region within which there is a 68% probability that a given PN belongs 
to the stellar population of M32, assuming a Gaussian stellar line-of-sight 
velocity distribution. The Vpn values are plotted against M5007 of the PNe 
to check whether bright and faint PNe, possibly powered by central stars 
having evolved via distinct binary and single-star channels, have different 
kinematics. A Spearman rank coefficient of 0.6 suggests a correlation in this 
figure, although this would not be dramatically significant since there is still 
a 2% probability that the plotted quantities follow each other monotonically 
only by mere chance. 



to ~ 10 3 - 10 4 yrs for stars leaving the HB red end) brown et al.l 
(2000) conclude that only a small fraction of HB stars in M32 live 
at its blue tip, so that in fact many more bright pAGB should be ob- 
served. The recent analysis of the UV colour-magnitude diagram 
(CMD) of M32 bv lBrown etal] l2008h further reinforces this re- 
sult. The findings of iBrown et alj ( I2000L l2008h not only contrast 
again with the notion that only hot and less massive HB stars should 
now form in the old stellar population of M32 (unless one consid- 
ers binary coalescence), but also add an additional problem for the 
current models of the late evolution of stars. 

To explain such a dea rth of pAGB stars, among other pos- 
sibilities |Bjownet_alJ (2000) suggested that pAGB could be ob- 
scured by circumstellar material, in particular as they cross the UV 
CMD before reaching the region where they peak in temperature 
at around 60,000 K. This is an interesting possibility considering 
that PNe ar e themselves affected by dust extinc tion and given that, 
assuming a Cardelli, C layton. & Mathisl dl989h reddening law and 
Rv = 3.1, the typical extinction of 0.7 magnitude at 5007 A for 
PNe in both star-forming and old galax ies ( Herrmann & Cia rdulld 



l2009tlRicher." Stasinka, & McCall 1999; |jacobv & Ciard ullo 1999) 
would translate into considerable reddening values of ~1.71 and 
~ 1 .67 in the HST far- and near-UV passbands (FUV and NUV), re- 
spectively. To actually check the impact of reddening in the pAGB 
population of M32, or at least for the fracti on that is presently 
powering a PN, we looked in the catalogue o f IBrown et al] feOOa. 
kindly provided by T. Brown) for the best candidate central star for 
our PNe sources and plotted their position in the UV CMD shown 
in Fig. [TO] The field-of-view of the lBrown et al.l HST observations 
falls entirely within the area covered by our SAURON data, being 
24" x 24" and only slightly offset from the center of M32, but un- 
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Figure 10. UV co lour-magnitude diagram (CMD) for the HST data of 
IBrown et alj |2008) with solid lines tracing, from top to bottom, the evo- 
lutionary path of pAGB, early-pAGB, and AGB-manque stars starting from 
the horizontal branch (shown by the grey thick line) and ending in the re- 
gion occupied by white dwarves. The asymptotic giant branch, from which 
pAGB and early-pAGB stars cross back the UV CMD is to the left and out- 
side the plotting range. The large circles indicate the location in this CMD 
of the brightes t and closes t UV sources to the PNe falling within the field- 
of-view of the lBrown et alj observations. Such PN central star candidates 
are colour-coded (from dark blue to white) ac cording to the brightness of 
the PN they could be powering. Assuming a Cardelli . Clayton. & Mathisl 
<1989h reddening law, the arrow at the top right of the figure shows the im- 
pact of dust reddening in the far- and near-UV HST passbands for an extinc- 
tion of 0.7 magnitude at 5007A. This reddening value for PNe is typically 
observed in both star-forming and old galaxies, and would suffice to bring 
the candidate central stars of the PNe in the optical regions of M32 close to 
the location where pAGB and early-pAGB stars cross the UV CMD before 
reaching the WD cooling curve, to the left of the bluer end of the horizontal 
branch. 



fortunately only 8 of our PNe end up within it. As we are looking 
to assess as conservatively as possible the impact of reddening on 
pAGB stars, in Fig. QJJ] we have located not only the closest but 
also the brightest of the UV stars near the location of our PNe. 
Except for PN 9, we could always find a pretty unique UV star 
within a radius of Of! 97 (half a FWHM) with a far-UV magnitude 
that would stick out from the distribution of magnitude values of all 
stars within 2" 5 (3 <tpsf) of our PN by at least 3 times the scatter 
in such values. Furthermore, except for PNe 12 and 15, such bright 
UV stars are always quite close to our PNe, within 0'.'25. For PN 
15 the relative large distance of the candidate star (0"63) is less of 
a concern, given that admittedly our Gaussian model does not fully 
reproduce the [O III] flux distribution of this source, but for PN 12 
(0"85) it is possible that the central PN star lies elsewhere near the 
PN and therefore also in the CMD of Fig. QJJ This is also most 
certainly the case of PN 9, for which the plotted position in Fig.UOl 
should be considered as upper limit for the far-UV brightness of its 
central star. It is important to note, however, that the central stars of 
even relatively faint PNe such as PN 9 and 12 should still be intrin- 
sically fairly bright pAGB or early-post AGB stars (epAGB) and, 
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in particular, that our PNe could not be powered by stars that have 
already rea ched the white-dwa rf (WD) cooling curve. Indeed, the 
models of Mari go et al ] d2004 see their Fig. 10) shows that at that 
stage it is possible to obtain only PNe of M5007 > 0.5, which are 
much fainter than the PNe that we can detect in M32. On the other 
hand, even considering the uncertain position of PN 9 and possibly 
PN 12 in Fig. [TO] we note that accounting for the typical reddening 
observed in PNe would bring back up and left the position of the 
central PN candidate stars of our brightest and faintest PNe near to 
the location in the UV CMD where pAGB and epAGB reach their 
peak temperature and spend most of their time, for FUV magni- 
tudes around ~ 20 and ~ 22, respectively, and a FUV-NUV colour 
~ —1.5. In fact, we note that PNe as bright as our PN 1, 5 and 
6, with A/5007 < —3.5, can only be powered by central stars that 
came from main-sequence stars of masses between 1.6 and 2.5 Mq 
and which would cross the UV CMD of Fig . 1 101 along the bright- 
est pAGB tracks (again, see Mar igo et al.l2004lT which implies that 
these sources actually must be reddened. 

We can double-check this result by using our integral-field 
data to estimate the position in the theoretical log T c g ~ log L 
Hertzsprung-Russell (HR) diagram of the PNe for which we can 
also detect the H/3 emission, and then compare their location to 
the evolutionary tracks of pAGB s tars. Through the models of 
Dopit a. Jacobv. & Vassiliadis] d 19921) it is indeed possible to esti- 
mate the temperature T e g and luminosity L of the PNe central star 
using the value for the [O IIl]/H/3 ratio and the H/3 flux, respec- 
tively. These models predict a conversion efficiency e between the 
bolometric luminosity L of the central star and the H/3 recombina- 
tion flux from the PN that is rather constant and independent of the 
gas temperature T a or metallicity Z gas . On the other hand, since 
Oxygen is the main coolant of PNe, the [O III] A5007 flux depends 
eavily on T c and Z gas and some knowledge of these quantities fur- 
ther helps constraining T e g with the [O IIl]/H/3 ratio. The restricted 
wavelength range of our SAURON data does not grant us access to 
the weak [O III]A4363 line that is normally used to measure T c , but 
we can use our spectroscopic data to measure the stellar metallicity 
along the line-of-sight of our PNe and thus gauge the metallicity of 
their central star, which in turn relates naturally to Z gas . 

There are 7 PNe in the central regions of M32 for which we 
could measure the [O IIl]/H/3 ratio (Tab. [JJ, with values ranging 
from to just above 6 to almost 17. The SAURON aperture spectra 
around these objects (extracted within a radius r — 3 <tpsf, Fig. [5} 
indicat e values for the H/3 abso rption line and the [MgFe50]' in- 
dex of lKuntschneretalJ feOld) between 1.77 - 1.92 A and 3.01 
- 3.08 A, respectively, which in turn yield an estimate for the 
stellar metallicity [Z/H] between -0.4 and -0.33 (see Fig. 3 of 
lKuntschneretalJl2010l). depending on whether one adopts the 
model s of ISchiavonl d2007h or of [Thomas, M araston. & Bendeil 
J2003h . With this gauge for the central star m etallicity we can infer 
from Fig. 3 of lDopita. Jacobv, & Vassil iadis a star temperature T e g 
around ~55,000 K for the PNe with low [Oni]/H/3 values ~ 6 - 
8, and above ~85,000 K for the PNe with [OlIl]/H/3> 11. From 
Fig. 1 of Dopita et al. this range of temperatures would then cor- 
respond to an average conversion efficiency e around 0.7%, which 
translate our observed H/3 fluxes into central star luminosities in 
the \og(L/ Lq) —3.0 - 3.2 range. These T e g and L values would 
place the 7 PNe with [O IIl]/H/3 measurements from our SAURON 
observations right below and to the right of the post-AGB evolu- 
tionary paths that are normaly used in the HR diagram to trace the 
luminosity evolution of PNe (see, e.g., Fig. 1 of Mendez et al. 2008 
and Fig. 3 of Schonberner et al. 2010), too cold to be on the cooling 
track of white dwarves (which cannot power PNe) and too faint to 



be on par with the pAGB stars that normally produce PNe (which 
range in luminosity between log(L/LQ) =3.5 - 4.0). Assuming 
a Solar or nearly Solar metallicity for the central stellar popula- 
tion of M32, as found bv lRose et al. (2005), would lead to an even 
more stringent result since in this case the central star T e ff of the 
highly-excited PNe would be better constrained to lower values by 
the models of Dopita et al. Thus, also for this second subsample of 
PNe it would seem necessary to invoke dust reddening in order to 
to reconcile our data with the model expectations, further support- 
ing the possible role of circumstellar dust in producing the apparent 
dearth of pAGB stars in M32. 



6 CONCLUSIONS AND FUTURE PROSPECTS 

Using SAURON data for the nearby elliptical M32 we have shown 
how, by means of a careful subtraction of the stellar background, 
integral-field spectroscopy (IFS) allows to detect the emission of 
single Planetary Nebulae (PNe) in the optical regions of early-type 
galaxies down to flux levels otherwise hardly accessible to standard 
narrow-band photometry. In turn, this makes it possible to trace the 
PNe luminosity function (PNLF) within a wider completeness limit 
and in galactic regions where stellar populations gradients can be 
the most extreme, and therefore where studying how the shape and 
normalization of the PNLF originate from the bulk of a galactic 
stellar population could be the most instructive. Moreover, we have 
shown that the incompleteness of a survey for PNe in the optical re- 
gion of a galaxy can be very well understood when using IFS data, 
and that the possibility to extract also the stellar kinematics further 
allows to check the membership of PNe on a local basis, that is, 
along the line-of-sight towards a given PN. Such a local comparison 
between the PNe and stellar kinematics could be useful to identify 
subpopulations of PNe of different origins. Finally, although with 
the SAURON data we could only confirm the high-ionisation nature 
of nearly half of our sources, on a longer wavelength range IFS has 
a great potential to study the nebular spectrum of extragalactic PN e, 
as already illustrated in the pioneering work of lRoth et al.l (12004). 

In the specific case of M32, with just two 10-minutes-long 
SAURON pointings, we have doubled the number of known PNe 
within the effective radius of this galaxy and detected PNe five 
times fainter than pre viously found with narrow-band imaging by 
ICiardulloetaIl fl989). while collecting slightly less the same num- 
ber of photons. Furthermore, accounting for the incompleteness of 
our survey across the entire luminosity range spanned by our de- 
tected PNe, we have concluded that the central PNe population of 
M32 is consistent with the generally adopted shape for the PNLF 
(that of Ciardullo et al.) and its typical normalization observed in 
early-type galaxies (i.e. with a PNe number density a of a few 10 -7 
per L©). Finally, thanks to exquisite HST data and to the proximity 
of M32 we were able to combine our PNe measurements with im- 
ages for the resolved UV-bright stellar population of M32 and con- 
clude that the PNe central stars may be considerably reddened by 
dust, a result that we could double-check thanks to the IFS nature 
of our data by estimating the temperature and luminosity for the 
central star of the PNe with [O Hl ]/H/3 measure ments. This finding 
would support the suggestion by iBrown et alj that circumnstellar 
dust extinction could explain the apparent dearth of pAGB stars in 
M32, a problem that may also apply to the halo of the Milky Way 
dWeston. Napiwotzki. & Catalanll201fj|) . 

Given the intrinsically small number of PNe that is expected 
in the central regions of such a small galaxy as M32, it is unlikely 
that with the current instruments (e.g., SAURON or VIMOS) fur- 
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ther IFS observations will reveal a sufficiently larger number of 
objects to add either to what found in this work or to what is al- 
ready known in the outskirts of M32 from narrow-band imaging or 
slit-less spectroscopy, even though more IFS data could allow for 
more emission-line diagnostics. 

On the other hand, there is a wealth of SAURON data for many 
more early- type galaxies that w as aquired over t he course of the 
SAUR ON ddeZeeuwet al]|2002h and ATLAS 3D JCappellari et al.l 
1201 lh surveys where the presence of PNe have not yet been sys- 
tematically investiga t ed. In fact, even though a few PNe were de- 
tected in ISarzi et all JioOrj) . many more PNe have certainly been 
missed because a large fraction of the data were spatially binned 
in order to extract a reliable stellar kinematics and due to the com- 
mon presence of diffuse ionised gas. Indeed spatial binning washes 
out the signal of the weaker PNe against the stellar background, 
whereas extended nebular emission makes it difficult to detect PNe 
unless their distinctive spectrum characterised by high-ionisation 
and narrow lines is purposely looked for. Paying attention also to 
the unresolved nature of the PNe emission and working with un- 
binned data using an educated guess for the stellar kinematics (for 
instanc e based on relati vely simple but functional dynamical mod- 
els, see lCap pellari 2008) it may be possible to tap a significant PNe 
population in the central regions of nearby early-type galaxies. In- 
deed, although the SAURON and ATLAS 3D surveys have targeted 
much more distant objects than M32, the ability to carefully sub- 
tract the stellar background will still enable the detection of fainter 
PNe in the central regions of these galaxies compared to narrow- 
band or slit-less spectroscopic surveys. Furthermore, the relatively 
large number of objects surveyed by these campaigns (260) will al- 
low to combine the measurements from different galaxies, and for 
instance to better explore the bright end of the PNLF or variations 
of its shape and normalization as a function of the stellar age and 
metallicity that will also be provided by the these SAURON data. 

To conclude, given the number of findings obtained here with 
a modest time investment (not even 30 minutes with overheads) on 
a 4m telescope, it is exciting to consider what could be learned 
on the PNe population in the optical regions of other galaxies, 
when the next generation of integral-field spectrographs will be 
mounted on 8m-c lass telescopes, such as MUSE on VLT (see, e.g., 
Bacon et al. 2006). For instance, with a longer wavelength range 
and a larger light bucket, MUSE will not only find many more and 
fainter PNe, but it will also allow full diagnostic of their nebular 
spectrum that could lead to a direct measurement of the gas metal- 
licity and to recognise other nebular so urces that could pose as a 
PNe in narrow-band surveys (see, e.g., iFrew & Parker 2010|). In 
turn the PNe metallicity could be compared with that of the un- 
derlying stellar population, which would also be better constrained 
thanks to an extended spectral range. Furthermore, the spatial res- 
olution that will be provided to MUSE by the use of adaptive op- 
tics will allow further comparisons with HST or other space-based 
measurements for the resolved stellar population of nearby galax- 
ies, such as those presented here for M32, which are bound to shed 
more light on the link between PNe and their parent stellar popula- 
tion. 
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